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There is a growing interest in very long baseline neutrino oscillation experimentation using ac- 
celerator produced neutrino beam as a machinery to probe the last three unmeasured neutrino 
oscillation parameters: the mixing angle 813, the possible CP violating phase 5c p and the mass 
hierarchy, namely, the sign of Am| 2 - Water Cherenkov detectors such as 1MB, Kamiokande and 
Super-Kamiokande have shown to be very successful at detecting neutrino interactions. Scaling up 
this technology may continue to provide the required performance for the next generation of exper- 
iments. This report presents the latest effort to demonstrate that a next generation (> 100 kton) 
water Cherenkov detector can be used effectively for the rather difficult task of detecting v e s from 
the neutrino oscillation —¥ v e despite the large expected potential background resulting from 7r°s 
produced via neutral current interactions. 
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I. INTRODUCTION 

It has been shown that there are certain advantages in 
a wideband neutrino beam over a now traditional narrow- 
band neutrino beam for neutrino oscillation experiments, 
especially for the observation of — > v e oscillation, pro- 
vided that the baseline is reasonably long (over 1,000 
km) 0. 

By broadening the range of energies with which multi- 
ple oscillations can be resolved, a richer parameter space 
can become accessible than is available to an experiment 
focused on only the first oscillation maximum. 

In the original paper proposing the use of a wide-band 
neutrino beam for a very long baseline neutrino oscil- 
lation (VLBNO) experiment [l] by the Brookhaven Na- 
tional Laboratory (BNL) neutrino working group (NWG) 
some simplifying assumptions were used. The sensitivi- 
ties presented in the paper relied on calculations based 
on 4-vector level Monte Carlo (MC) simulations, a sim- 
ple model for energy resolution and certain assumptions 
about reconstruction capabilities. Namely, a detailed 
detector simulation for the proposed water Cherenkov 
detector was not included. In addition it was assumed 
that the signal events were only from quasi-elastic (QE) 
charged current (CC) scattering (v e + n — > e~ + p) and 
the background events were only from single 7r° neutral 
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current (NC) interactions v + N^v + n Q + N' where N 
and N' are nucleons. 

In order to gain a better insight on this idea of the 
VLBNO experiment with a wide-band neutrino beam for 
— > v e oscillation, we performed a more sophisticated 
and elaborate multivariate likelihood based analysis us- 
ing a full MC simulation that included inelastic neu- 
trino interactions, water Cherenkov detector response, 
and well-tuned event reconstruction algorithms. This 
MC simulation and reconstruction programs were devel- 
oped and fine-tuned for the Super-Kamiokande-I experi- 
ment (SK-I) 0. 



II. ENERGY RECONSTRUCTION AND EVENT 
DESCRIPTION 

For the discussions in this paper, we calculate a recon- 
structed neutrino energy using the formula 

E _ m N E e 

rec rriN — (1 — cos9 e )E e ' 

where my, E e and e arc the nucleon mass, the recoil 
electron energy and the scattering angle of the recoil elec- 
tron with respect to the incident neutrino beam, respec- 
tively. In a strict sense, this quantity represents the in- 
cident neutrino energy only when the event is produced 
by CC QE scattering and the Fermi motion of target 
nucleons is ignored. Nonetheless, we show in Figure [1] 
the neutrino energy and E rec for CC QE (top) and all 
CC (bottom) events to compare the two energy distribu- 
tions for the two classes of CC events. Although E rec for 
the CC events does not reproduce the incident neutrino 
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FIG. 1: The distributions of neutrino energy (solid lines) and 
reconstructed neutrino energy (dotted lines) of single ring e- 
like events originating from CC QE interactions (top figure) 
and from all CC interactions (bottom figure). Note that a 
dip at about 2 GeV is due to — > v e neutrino oscillation. 
The open square boxes with error bars indicate statistical 
uncertainties. 



spectrum as well as it does for CC QE events, it still 
reproduces it quite well. 

The initial event selection is made to maximize the 
number of interactions that are consistent with coming 
from v e appearing at the far end of the baseline in the 
beam while minimizing events such as CC or any 
events from NC interactions. Needless to say, the clear- 
est event signature comes from v e CC QE interactions. 
Inelastic v e CC events also carry information about os- 



cillations but are more difficult to cleanly select and have 
somewhat worse energy resolution as shown in Figure [T] 
Regardless of the interaction, a recoil proton will rarely 
be above its Cherenkov threshold of about 1.25 GeV/c 
and will not contribute notably to the event signature. 
For these reasons, a v e appearance candidate event will 
be initially selected if its signature is consistent with be- 
ing a Cherenkov radiation pattern from a single electron 
(termed single ring, e-like). 

A CC QE event will usually be distinguishable from 
one arising from a v e CC QE interaction as the result- 
ing Cherenkov ring produced by will have a sharper 
rise and fall in the amount of Cherenkov light at the 
edges of the ring than one due to an electron. This is 
because a muon presents a relatively straight line source 
of Cherenkov light while an electron initiates an elec- 
tromagnetic shower. This shower consists of a distribu- 
tion of track directions due to the number of particles in 
the shower and that they experience larger multiple scat- 
tering. As the muon energy approaches the Cherenkov 
threshold there is an increasing chance that multiple 
scattering and the collapsing Cherenkov light cone will 
smooth the rise of the edge of the ring and potentially 
mimic an electron. 

Likewise, inelastic and NC interactions with a non-zero 
pion multiplicity can have additional light that is not con- 
sistent with the pattern of a single e-like ring. Depending 
on the exact event topology these events are either easily 
ruled out or present challenging backgrounds. Charged 
pions above their Cherenkov threshold of 160 MeV/c pro- 
duce a /i-like ring and are rejected. Events with visible 
light from both pions and the primary charged lepton are 
strongly rejected based on ring counting. 

Events with a final state consisting of only a single neu- 
tral pion present a particular challenge. The ir° decays to 
a pair of gamma rays, each of which initiate an electro- 
magnetic shower and produce an e-like ring of varying in- 
tensity and direction. These ir° events can become back- 
ground depending on how the 7r° decays to two gamma 
rays. At one extreme, the decay is symmetric such that 
the gamma rays have similar energies. If the original ir° is 
boosted enough, the two gamma rays are nearly collinear 
and produce mostly overlapping rings. They can then be 
impossible to be distinguished single electron events with 
an energy equal to the sum of the two gammas. At the 
other extreme, a highly asymmetric ir° decay results in 
one strongly boosted and one strongly retarded gamma 
ray. In the lab frame, the higher energy gamma will pro- 
duce a single e-like ring while the other may produce no 
discernible light. 

The events with a ir° in the final state can come from 
two sources: (1) NC interactions and (2) charge ex- 
changes of charged pions inside the target nucleus or with 
an oxygen or a hydrogen nucleus while traveling in water. 

Finally, there are events from electron neutrino inter- 
actions that are not from v e that appeared in the beam 
from neutrino oscillation v„ — > v e . Rather, there is an 
intrinsic v e component in the unoscillated neutrino beam 
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originating from muon or kaon decays. The beam used in 
this study has an intrinsic v e component which is 0.7% of 
the flux. The background events from these neutrinos 
are irreducible. 

In this analysis, to classify the events for the initial se- 
lection, the reconstruction algorithms used for the SK-I 
atmospheric neutrino analysis are employed. In addition, 
a special algorithm called POLfit [8] is used. This algo- 
rithm has been found to be extremely useful for removing 
single 7T° background events in samples containing events 
that are classified as single-ring and e-like @ by the stan- 
dard SK-1 reconstruction. This paper describes how in- 
formation provided by POLfit, together with other useful 
variables, can be effectively used to significantly reduce 
the 7r° background while retaining a sufficient efficiency 
for the signal. 



III. MONTE CARLO EVENT SAMPLE 

The Monte Carlo event sample used in this study was 
originally produced by the SK-I collaboration to simulate 
atmospheric neutrino events detected by the SK-I detec- 
tor and corresponds to about 100 years of the exposure. 
The detailed description of the Monte Carlo generation as 
well as the event reconstruction is described elsewhere Q . 

The energy spectra for za, and v e produced in the 
Earth's atmosphere are different from those in the neu- 
trino beam proposed by the BNL NWG. To account for 
this, each event is given a weight based on the neu- 
trino energy such that the shape of the resulting spec- 
trum matches that expected from the beam. An addi- 
tional normalization is applied so that the total number 
of CC QE events is as expected. 

Assuming no oscillations, it is expected that there will 
be 12,000 such events for a detector with SK-I efficien- 
cies and fiducial volume of 500 kton (22.2 x SK-I) which 
is placed 2,540 km from the neutrino production target 
and for five years of running as described in the paper by 
BNL NWG When oscillations are applied, this sam- 
ple is weighted by the oscillation probability taking into 
account full three-neutrino mixing and matter effects. 

The oscillated v e sample is prepared by similarly 
weighting and normalizing the v e atmospheric MC events 
to what is expected from the component of the beam, 
that is, under the assumption of 100% — » v e oscilla- 
tion. Actual oscillation parameters are then applied to 
this sample by weighting to the appropriate oscillation 
probability. The intrinsic v e component in the beam is 
taken to be the shape predicted by the beam MC simu- 
lation and is normalized to be 0.7% that of the com- 
ponent. To include the effect of neutrino oscillation, the 
contribution of the v e component is calculated only when 
v e survives as v e at the water Cherenkov detector placed 
at a given baseline. 

The expected numbers of the events to be observed 
from the intrinsic anti-z/^ component of the beam are es- 
timated to be 2.2% and 2.5% of the expected numbers 



of the background events from the (background- 1) for 
the baselines of 2,540 km and of 1,480 km, respectively 
(see the next section for the description about the base- 
line) . The expected numbers of the events to be observed 
from anti-z^ e component of the beam are estimated to be 
less than 4.5% and 4.3% of the expected numbers of the 
background events from the intrinsic v e component of 
the beam (background-2) for the baselines of 2,540 km 
and of 1,480 km, respectively. Since these contributions 
are much smaller than the events from the background- 1 
and the background-2, although not negligible, we will 
not take into account them further in this report. 

To select MC v e appearance candidate events an initial 
set of cuts based on the standard SK-I codes is used: (a) 
there is one and only one e-like ring, (b) the reconstructed 
event vertex is at least 2 m from any inner photomulti- 
plier (PMT) surface and (c) there is no activity in the 
outer (veto) detector. 

In addition to these standard cuts, events produced by 
neutrinos with energies greater than 10 GeV are ignored 
due to lack of statistics in the atmospheric MC sample. 
The beam has a tail that extends up to 15 GeV but its 
contribution to the background is negligible as will be 
shown below. 



IV. NEUTRINO OSCILLATION 
PROBABILITIES 

The oscillation probabilities applied to the components 
of the beam neutrino flux use the following parameters: 

Am^ = 7.3 x 10- 5 eV 2 , 

Am§ x = 2.5 x 10~ 3 eU 2 , 
sin 2 26>i2 = 0.86, 
sin 2 2023 = 1.0, 
sin 2 26»i3 = 0.04, and 

S C p = 0°, ±45°, ±135°. 

Note that the sign convention of Sep follows that of 
the paper by BNL NWG Q. Unless otherwise stated, in 
this paper the value of Sep is +45° and the baseline is 
2,540 km which corresponds to the distance from BNL 
to Homestake Mine in South Dakota. Later we will de- 
scribe the results with different values for Sep as well 
as with the baseline of 1,480 km which corresponds to 
the distance from Fermilab to Henderson Mine in Col- 
orado. The result with the baseline of 1,480 km was for 
the UNO detector proposed first in 1999 and later to be 
built in Henderson Mine. See the references [H for the 
detail of the UNO detector. These results are similar to 
what is expected for the Fermilab to Homestake baseline 
of 1,300km. This shorter baseline is considered by apply- 
ing inverse-square scaling of the flux and recalculation of 
oscillation probabilities. 

The neutrino oscillation probabilities were calculated 
by the nuosc || package. The calculator assumes 3-z^ os- 
cillation with possible CP symmetry violation and mat- 
ter effects. The matter densities can be vacuum, non-zero 
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constant or stepwise-constant matter densities. The non- 
constant PREM [6| and arbitrary user provided density 
profiles are also supported. It calculates the probabili- 
ties either analytically for constant matter densities or 
for any density profile it can use a 5th order adaptive 
Runge-Kutta stepper. The package has been validated 
against published calculations 0). 



V. POLFIT 

In this section we briefly describe the w° fitter called 
Pattern Of Light fit {POLfit) POLfit is applied to 
events which are initially classified as single-ring and e- 
like by the standard SK-I analysis codes in order to iden- 
tify 7T° background events in this event sample. It as- 
sumes the event is a 7r n event and that the ring found 
by the standard codes is one of the gammas from the 
7r° decay. It will then determine the direction and en- 
ergy of a secondary gamma that, along with the primary 
one, is most consistent with the pattern of light collected 
by the PMTs. POLfit produces two possible secondary 
gammas by employing two algorithms. One is optimized 
assuming the secondary ring is overlapping with the pri- 
mary ring (forward-algorithm). The other is optimized 
to find a second ring in the wider angular region (wide- 
algorithm) . 

In each algorithm, the pattern of light expected from 
the full 7T° decay is calculated from a set of templates and 
compared to the collected light and a likelihood is formed. 
These templates are pre-determined using the detector 
simulation. The energy and direction phase-space of the 
secondary gamma is searched to find the point that leads 
to the maximum likelihood. Each algorithm supplies the 
direction of two photons, their energies, the two photon 
invariant mass, and the maximum likelihood value ob- 
tained. It was found that the wide- algorithm method 
more reliably finds a real second ring more often than 
the forward-algorithm. Therefore, in this report we use 
mostly the information provided by the wide-algorithm 
unless otherwise stated. 

To demonstrate the power of POLfit, Figure [2] shows 
the relative tt° reconstruction efficiencies with (solid cir- 
cles) and without (open circles) POLfit employed as a 
function of the opening angle between the two recon- 
structed gammas from the ir° decay in the lab frame and 
for single ir° events produced by the NC interactions. To 
be accepted in this sample, the event must have a two- 
photon invariant mass (m 77 ) within 2er of the expected 
value given the event came from a ir°. Here, the neutrino 
energy spectrum is that of the original SK-I atmospheric 
muon neutrinos. In the case without POLfit, only events 
identified as having two e-like rings by the standard SK-I 
software can be candidates to be a 7r°. With POLfit in- 
formation, 7T° events that are classified as single-ring by 
the standard selection may now be selected. 

Depletion of ir° detection efficiencies at a small open- 
ing angle is due to overlap of the two e-like rings. At 
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FIG. 2: 7T° reconstruction efficiency with the standard SK- 
I codes (open circles) and after accepting additional events 
with missing gamma found by POLfit (solid circles). 



large opening angle it is due to the second ring being too 
faint. As clearly seen in Figure [2] POLfit improves the 
7T° reconstruction efficiency significantly. 



VI. USEFUL VARIABLES TO DISTINGUISH 
THE SIGNAL FROM THE BACKGROUND 

A cut on m 77 reduces background from single ir° NC 
events but it also reduces signal efficiency somewhat. 
Alone, it is not enough to reduce the background to a 
level comparable to the oscillation signal. It is, there- 
fore, desirable to find additional distinguishing features. 

In this section, nine such features (variables) including 
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are described. Distributions of these variables are 
shown for events that are oscillated v e CC (labeled "sig- 
nal" and with solid lines) and NC, CC and beam-i/ e 
CC events ("background", dotted lines). The distribu- 
tions are plotted for different E rec regions in steps of 0.5 
GeV from to 2 GeV, 2 GcV < E rec < 3 GeV and E rec > 
3 GeV. 

In general, at lower energies (sub-GeV range), the 
background is largely from the misidentified single-7r° 
events discussed above. At higher energies, non-QE in- 
teractions become dominant leading to additional ways 
for background to mimic v e CC QE events. 



A. Reconstructed it mass (m 77 ) 

As shown above, this variable is useful to remove back- 
ground events that come from single ir° NC production. 
In Figure |3] the distributions of shown. In the 

lower reconstructed neutrino energy region (E rec < 2 
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FIG. 3: The distributions of reconstructed ty° mass of single e- 
like ring events for signal (solid line) and background (dotted 
line). 



FIG. 4: The distributions of the energy fraction of the second 
ring found by POLfit in 1-ring events for the signal (solid line) 
and background (dotted line). 



GeV), the distributions show a prominent ir° peak. At 
higher energies, the contribution from multi-pion produc- 
tion increases and POLfit begins to reconstruct the cor- 
rect second gamma poorly or returns an arbitrary result 
that is not associated with any physical particle. Because 
of this, the tt° peak disappears in the higher energy re- 
gion. 



B. Fraction of energy in the second ring (E/, ac ) 

In general, the second e-like ring found by POLfit has 
less energy than the one found by the standard SK-I re- 
construction software. Furthermore, POLfit must find a 
second ring and if the event does not in fact have one, 
this wrongly reconstructed ring tends to have much less 
energy than the primary one. This effect is shown in Fig- 
ure S] where the ratios of the energy of the second ring to 
the sum of two ring energies are plotted. 



C. Difference in log it likelihood ratio (Alog 7r°-lh) 

As mentioned above, POLfit employs two algorithms, 
the wide- and forward-algorithm. Each provides a log- 
likelihood that the event is a 7r°. Signal events actually 
contain only a single ring. Nevertheless, with such events, 
POLfit will return a second, fake ring. In such cases, both 
algorithms tend to place this fake ring in the vicinity of 
the single primary ring. Therefore, the two algorithms 
tend to find similar fake rings. This makes the two like- 
lihoods similar to each other and the log-likelihood ratio 



(difference between two log-likelihoods) tends to be sym- 
metric around zero. On the other hand, in the case of a 7r° 
background event the secondary photon from a 7r° decay 
is not necessarily confined in the direction of the primary 
photon. This makes the log-likelihood ratio asymmetric 
with a long left tail, especially at lower energies. The dif- 
ference is significant in the lower energy region as shown 
in Figure [5] However, this trend changes above E rec — 
1.0 GeV where the contribution from multi-pion produc- 
tion starts to increase and dilutes this effect. 



D. Direction cosine of the e-like ring (cos#) 

The direction cosine of the primary e-like ring with re- 
spect to the neutrino beam is a good discriminator to sep- 
arate the signal from the background. Its distributions 
are shown in Figure [6j This variable does not depend on 
the information from POLfit. 

As the energy of the neutrino decreases, the outgoing 
electron direction is distributed over a wider angle with 
respect to the neutrino direction. As the energy becomes 
comparable to that of Fermi motion, the distribution is 
further widened because the momentum of the target nu- 
cleon becomes significant and will tend to randomize the 
electrons direction. This latter effect is less for the more 
massive 7r°. 
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FIG. 5: The distributions of the difference in log n likelihood 
between the two algorithms. The larger this variable, the 
more likely an event is the 7r° background. The distributions 
in solid lines are for the signal and those in dotted lines are 
for the background. 



FIG. 6: The distributions of the directional cosine of the pri- 
mary e-like ring with respect to the neutrino beam direction. 
The distributions in solid lines are for the signal and that in 
dotted lines are for the background. 



E. Total charge to ring energy ratio (Q/E) 



Light collected in an event that is not consistent with 
the single ring found by the standard SK-I codes is an in- 
dication of unreconstructed or missed particles and thus 
a non-QE interaction. A measure of this light is the ra- 
tio of the total charge, in unit of PE (photoelectron), 
collected by all PMTs to the energy associated with the 
primary reconstructed ring in MeV. Some background 
events where only one ring is found are expected to pro- 
duce some light which is not identified as a ring. The 
distributions of this variable are shown in Figure [71 The 
separation is most pronounced at lower reconstructed en- 
ergy where any 7r°s are likely to be less boosted and not 
put all their light in the forward direction. Thus in an 
asymmetric ir° decay the second ring tends to be weaker 
than the standard SK codes can detect. 
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F. log particle-identity-likelihood (log pid-lh) 

The standard SK-I reconstruction software provides 
likelihoods of a Cherenkov ring to be e-like (L e ) and fj,- 
like (Lfj). From these two likelihoods, the logarithm of 
the ratio of likelihoods, log(L M /L e ) (log pid-lh) is used as 
a good measure to separate electrons from muons. The 
more negative log pid-lh is, the more e-like a Cherenkov 
ring is. When two e-like rings, such as the photons from 
a 7T° decay, overlap and reconstruct as a single ring they 
can produce a log pid-lh that is more e-like than a single 



FIG. 7: The distributions of the ratio, the total charge in PE 
to the ring energy in MeV are shown. The distributions in 
solid line are for the signal and those in dotted line are for 
the background. 



electron at a comparable energy would have. Figure [8] 
shows this small but noticeable effect. 
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FIG. 8: The distributions of the log-likelihood ratio between 
e-like and /i-like of the primary e-like ring. The distributions 
in solid line are for the signal and those in dotted line are for 
the background. 



FIG. 9: The distributions of the log 7r°-likelihood of single e- 
like ring events. The distribution in solid line is for the signal 
and that in dotted line is for the background. 



G. log 7r°-likelihood (log 7r°-lh) 

POLfit provides the logarithm of a likelihood that the 
event consists of a single 7r°. For events that actually 
do consist of a single 7r° this variable is expected to be 
smaller (more negative), namely more 7r°-like, than that 
for a signal event. This trend shows up in lower energy 
region (E rec <1 GeV) as shown in Figure[5J In the higher 
energy region, the distribution tends to be narrower for 
the signal events. 



H. Cherenkov angle (Cangle) 

The distribution of reconstructed Cherenkov angle is 
expected to be different between the signal and back- 
ground events and it depends on whether there is overlap 
between two Cherenkov rings and on the energy of the 
primary ring. In Figure [10] the Cherenkov angle distri- 
butions for different E rec regions are shown. The shape 
of the distribution for the signal events differs from that 
for the background events in most of the energy regions, 
although degrees of differences vary from energy region 
to region. 



I. Log ring-count likelihood ratio (Alog ring-lh) 

To decide how many Cherenkov rings there are in an 
event, two hypotheses are compared: the first hypothesis 
is that there is only one Cherenkov ring and the second is 
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FIG. 10: The distributions of the measured Cherenkov angle 
of the primary e-like ring. The distribution in solid line is for 
the signal and that in dotted line is for the background. 



that there is an additional ring in the event. The compar- 
ison is made using log-likelihood ratio, calculated by the 
standard SK-I codes, of the two hypotheses. In Figure [TT1 
this log-likelihood ratio is plotted for different E rec re- 
gions. The shape of the distribution for the signal events 
differs from that for the background events in most of 
the energy regions, although degrees of differences again 
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FIG. 11: The distributions of log-likelihood ratio of two hy- 
potheses on the number of Cherenkov ring in an event: there 
is only one Cherenkov ring or there are more. If the log- 
likelihood ratio is negative, the event is less likely to have 
an additional ring. The distributions in solid line are for the 
signal and that in dotted line are for the background. 



vary over the different energy regions. 



VII. DISCRIMINATOR: LIKELIHOOD 
FUNCTION AND LIKELIHOOD RATIO 

A series of cuts on the variables described in the pre- 
vious section can reduce the background well. However, 
too many cuts in series reduce the signal selection effi- 
ciency. The problem for this analysis is that none of the 
variables alone can separate the signal from the back- 
ground powerfully. Thus, this situation is well suited 
for a multivariate analysis where a likelihood function 
is defined and utilizes all variables that have noticeable 
distinguishing power. We have shown that the nine vari- 
ables described above are distributed differently depend- 
ing on the source of the events (signal or background), 
although the differences may not be large. As we will 
see, an accumulation of relatively small differences, when 
the correlations among the variables are small, can make 
a bigger difference. 

From each distribution shown in Figures [51 1111 a 
probability is calculated for an event to have a value 
of the corresponding variable i as signal pf and as 
background p\. Then the log- likelihood is defined by 
log(Z s ) = Ei=i,..,9 log(p?) as signal and by \og{L b ) = 
£i=i,..,9 log(Pi) as background. These probabilities are 
calculated in the same bins of E rec as used in Figures [3]- 

Then the difference between two log-likelihoods, 



FIG. 12: The distributions of log-likelihood ratio of two hy- 
potheses on the origin of events: signal vs. background. The 
distributions in solid line are for the signal and those in dotted 
line are for the background. 

Alog(L) = log(ib) — \og{L s ), is calculated to decide 
whether to accept or reject the event. Figure [T^] shows 
A log(L) distributions for the different E rec regions for 
the signal events (solid line) and the background events 
(dotted line). The smaller Alog(L) is, the more likely an 
event is a signal event. It is clearly seen that the A \og(L) 
distribution of the signal events differ significantly from 
that of background events over wide range of E rec . 

VIII. RECONSTRUCTED NEUTRINO ENERGY 
DISTRIBUTIONS WITH A CUT ON Alog(L) 

The final acceptance of an event is determined by its 
value of Alog(L). As the distributions of the nine vari- 
ables used to define the likelihood depend on E rec , so 
does the distribution of A log(L). Therefore, in order not 
to change the energy spectrum unnecessarily, we adopt a 
strategy to keep the signal detection efficiency constant 
over a wide range of E rec by changing the A log(L) cut 
according to E rec . The cut is selected to retain 40% of 
the signal that passes the standard SK-I cuts for all E rec . 

In presenting the results, first we illustrate the per- 
formance of the traditional analysis represented by the 
standard SK-I codes for comparison. Figure [13] (top) 
shows the E rec distributions of the signal (dashed line), 
of the background (background- 1) mostly from NC inter- 
actions (dotted line), and of the irreducible background 
(background-2) from the v e contamination in the neu- 
trino beam (dash-dotted line). The signal is overwhelmed 
by the background, especially in the low energy region. 
In this case, we find 700 signal events, 1,877 background 



All event 
Signal 

Background from n c and other events 
Background from v t beam 




TABLE I: Summary of the numbers of the signal and back- 
ground events, and the signal-to-background ratio using 
events from background- 1 with the cut on Alog(L) to retain 
40% of the signal events. The signal-to-background ratio with 
one variable removed is compared to the ratio with all nine 
variables. Note that the numbers of the signal (background- 
2) events are not always 280 (45) due to the finite bin size 
used for the distribution of the variable removed. 



Variable removed 


Sig 


Bkg-1 


Bkg-2 


Sig/Bkg-1 


None 


280 


87 


45 


3.22 


A log 7T U -lh 


281 


102 


45 


2.75 


Q/E 


281 


94 


45 


2.98 


log 7T°-lh 


278 


94 


47 


2.98 


log pid-lh 


277 


94 


42 


2.96 


E^ rac 


281 


98 


45 


2.85 




280 


105 


45 


2.66 


cosf? 


279 


101 


45 


2.76 


C angle 


280 


98 


45 


2.86 


A log ring-lh 


277 


95 


45 


2.93 



All events 
Signal 

Background from tt c and other events 
Background from V, beam 
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FIG. 13: Top: The distributions of the reconstructed neutrino 
energy with the standard Super-Kamiokande cuts. The dis- 
tributions in dashed line are for the signal and those in dotted 
(dash-dotted) line are for the background-l(-2). Sep = +45° 
and the baseline is 2,540 km. Bottom: The same distributions 
except that, in addition to the standard Super-Kamiokande 
cuts, the cut on AlogL is applied. 



events from background- 1, and 127 background events 
from background-2. 

On the contrary, if we retain 40% of the signal events 
after the cut on A log(L) introduced in this paper, the 
background- 1 events are strongly suppressed as shown in 
Figure [131 (bottom) . We find 280 signal, 87 background- 
1, and 45 background-2 events. 

The analysis presented here uses nine variables to 
define the likelihood. The usefulness of all nine vari- 



TABLE II: Summary of the numbers of the signal and back- 
ground events for different values of Scp using events from 
background-1 (Bkg-1) and background-2 (Bkg-2) with the cut 
on A log(L) chosen to retain 40% of the signal events. The er- 
rors are statistical due to the limited size of the Monte Carlo 
event sample. 



Scp 


Sig 


Bkg-1 


Bkg-2 


+135° 


386+6 


89+8 


45+1 


+45° 


280+5 


87+8 


44+1 


0° 


197+4 


90+8 


44+1 


-45° 


159+3 


87+8 


44+1 


-135° 


263+3 


87+8 


45+1 



ables is evaluated by examining how much the signal-to- 
background (S/B) ratio changes when each variable is 
removed in turn from the overall likelihood function. Ta- 
ble [I] summarizes the numbers of events and the S/B ratio 
for each case where only the contribution to background- 
1 is considered. The S/B ratio is always appreciably 
smaller when one of the nine variables is removed from 
the likelihood function than when all nine are included. 
Each variable helps to enhance the signal. 



IX. E rec DISTRIBUTIONS AND CP- VIOLATING 
PHASE S C p 

In the previous section we showed that with a set of 
appropriate cuts, the background contribution can be 
suppressed down to a reasonable level, while retaining 
enough statistics for the signal. In this section, we see 
whether this set of the cuts is still useful for other values 
of the CP-violating phase Sep- Table HIl lists the numbers 
of events from the signal, background-1 and background- 
2 for various values of Scp- 
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FIG. 14: The distributions of the reconstructed neutrino en- 
ergy is shown for events that pass the standard SK-I cuts and 
the cut on A log(L) such that 40% of the signal events are re- 
tained. The distributions in dashed line are for the signal and 
those in dotted (dash-dotted) line are for the background- 1 
(background-2). Sep = +45° and the baseline is 1,480 km. 



TABLE III: Summary of the numbers of signal and back- 
ground events for different values for Sep using events from 
background-1 (Bkg-1) and background-2 (Bkg-2) with the cut 
on Alog(L) to retain 40% of the signal events and using a 
baseline of 1,480 km. The errors are statistical due to the 
limited size of the Monte Carlo event sample. 



Scp 


Sig 


Bkg-1 


Bkg-2 


+135° 


646±9 


238±23 


133±2 


+45° 


698±11 


235±23 


133±2 


0° 


498±8 


230±23 


132±2 


-45° 


356±6 


250±24 


134±2 


-135° 


609±9 


237±23 


134±2 



X. E rec DISTRIBUTIONS AND BASELINE 



It is interesting to see how the baseline will change 
the results of similar analyses presented in the pre- 
ceding section. For this study, the cut on Alog(L) 
is used again to retain 40% of signal and the base- 
line is assumed to be 1,480 km (Fermilab to Hender- 
son Mine). Figure [H] shows the E rec distribution of 
the signal (dashed line), background-1 (dotted line) and 
background-2 (dash-dotted) for 5 = +45°. Table Hill lists 
the numbers of events from the signal, background-1 and 
background-2 for different values of 5c p- 



XL SOURCES OF BACKGROUND EVENTS 

It is important to know where the background events 
come from. Information such as the true energy of neutri- 
nos and nature of interaction of neutrinos that produce 
the background events is very useful for design of the 
neutrino beam for a VLBNO experiment. 

Figures [T5JTJJ] show the neutrino energy distributions of 
the signal, background-1 and background-2 events for the 
baseline of 2,540 km and 1,480 km, respectively. These 
events are chosen with the A log(L) cut that retains 40% 
of the signal events after the first set of the cuts (the 
standard SK-I cuts). 

As mentioned earlier, for the most of the results pre- 
sented in this report, we apply the cut on the neutrino 
energy at 10 GeV. To justify this cut, we checked to see 
how much more the background contribution would in- 
crease if we allowed events produced by neutrinos whose 
energies were greater than 10 GeV and less than 15 GeV. 
For these additional events we fixed the weight value to 
that used at 10 GeV. The atmospheric v e flux times the 
v e cross section decreased by about 60% while that for 
the wideband beam by about 50% in this energy re- 
gion. Therefore the atmospheric v e and spectrum are 
only slightly softer than the wideband spectrum from 
10 GeV to 15 GeV. This argument, therefore, gives a 
reasonable estimate of the extra contribution from high 
energy neutrinos. All the neutrino oscillation parameters 
are the same as in the case with Sep — +45°, the base- 
line of 2,450 km, and the 40% A log(L) cut. The numbers 
of the signal, background-1, and background-2 events ac- 
cepted are 281, 91 and 45, respectively, which should be 
compared with 280 for the signal, 87 for the background- 
1 and 45 for the background-2 with the neutrino energy 
cut at 10 GeV. Thus the percent increases in the number 
of the signal, background-1 and background-2 events are 
1%, 4%, and 0%, respectively. 

Tables IIVI and [V] summarize all neutrino interactions 
that produce candidate events passing all cuts with 40% 
efficiency for the baseline of 2,540 km and 1,480 km, re- 
spectively. In these tables, tt ,tt ± , and nw stand for sin- 
gle 7T°, single 7r ± and multiple n production, respectively. 
The C07T label stands for tt production via coherent inter- 
action with the oxygen nucleus. The "others" includes 
eta and kaon production as well as deep inelastic scatter- 
ing (DIS). However, for the background, the label "nir" 
includes DIS. Note that in the tables there is no contribu- 
tion From NC interaction for signal events by definition. 
CC QE contributions to background-1 are due to in- 
teractions. Since the cross section of the coherent pion 
production is known only to an accuracy of ±30% (see 
reference 3), we estimate the effect of this error on our 
background-1 contribution by varying the cross section 
by ±30%. For S C p = ±45° with the baseline of 2,450 km 
±30% changes in this cross section result in changes in 
the numbers of background-1 events by ±3%. 



11 




E, (MeV) 



FIG. 15: The distribution of the energies of neutrinos that 
produce the signal, background- 1 and background-2 events. 
In addition to the standard Super-Kamiokande cuts, the cut 
on Alog(i) is applied in such a way to retain 40% of the 
signal events that survive by this cut. The distribution in 
dashed line is for the signal and that in dotted (dash-dotted) 
is for the background-l(-2). Sqp — +45° and the baseline is 
2,540 km. 
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TABLE IV: The percent contributions from events pro- 
duced by different interactions for signal events and for the 
background- 1 events are summarized for Sep = +45° with 
the cut on A log(L) to retain 40% of the signal events and 
with the baseline of 2,540 km. 



Interaction 


Erec range (GeV) 


Sig 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-3.0 


3.0- 


CC QE 

CC 7T° 
CC TT* 

CC n-K 
CC others 


86% 
2% 

11% 
1% 
07o 


79% 
3% 
15% 

3% 

noy 
07o 


63% 
4% 

28% 
3% 

1 07 
170 


82% 
2% 
13% 

4% 

noy 
O/o 


28% 
5% 
37% 
25% 

907 

3 7o 


46% 
6% 
30% 
14% 

O07 


Bkg-1 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-3.0 


3.0- 


CC QE 

CC K° 

CC 7r± 
CC mr 
CC others 

NC 7T° 
NC TT* 

NC nK 

NC C07T° 

NC elastic 


7% 
0% 
2% 
0% 
0% 
23% 
64% 
0% 
5% 
0% 


5% 
1% 
5% 
0% 
0% 
53% 
10% 
13% 
14% 
0% 


6% 
4% 
0% 
3% 
0% 

60% 
6% 
5% 

16% 
1% 


0% 
6% 
1% 
6% 
0% 
59% 
0% 
16% 
12% 
0% 


0% 
4% 
0% 
0% 
21% 
18% 
0% 
55% 
2% 
0% 


0% 
0% 
0% 
0% 
3% 
0% 
0% 
92% 
5% 
0% 



TABLE V: The percent contributions from events produced 
by different interactions for the signal events and for the 
background- 1 events are summarized for Sep — +45° with 
the cut on A log(L) to retain 40% of the signal events and 
with the baseline of 1,480 km. 



Interaction 




Erec rang 


e (GeV) 






Sig 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-3.0 


3.0- 


CC QE 


82% 


77% 


52% 


42% 


54% 


53% 


CC 7T° 


3% 


3% 


6% 


7% 


4% 


5% 


CC 7r± 


14% 


18% 


30% 


36% 


30% 


30% 


C nK 


1% 


2% 


12% 


13% 


11% 


10% 


CC others 


0% 


0% 


0% 


2% 


1% 


2% 


Bkg-1 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-3.0 


3.0- 


CC QE 


8% 


4% 


1% 


3% 


0% 


0% 


CC K° 


0% 


2% 


6% 


4% 


0% 


0% 


CC 7r± 


6% 


6% 


2% 


1% 


1% 


0% 


CC nK 


0% 


0% 


1% 


1% 


0% 


0% 


CC others 


0% 


0% 


0% 


4% 


2% 


7% 


NC 7T° 


23% 


58% 


61% 


29% 


23% 


0% 


NC TT* 


59% 


5% 


6% 


0% 


0% 


0% 


NC nK 


0% 


11% 


5% 


23% 


65% 


93% 


NC C07T° 


4% 


14% 


17% 


35% 


9% 


0% 


NC elastic 


0% 


0% 


1% 


0% 


0% 


0% 



FIG. 16: The distribution of the energies of neutrinos that 
produced the signal, background- 1 and background-2 events. 
In addition to the standard Super-Kamiokande cuts, the cut 
on Alog(L) is applied in such a way to retain 40% of the 
signal events that survive by this cut. The distribution in 
dashed line is for the signal and that in dotted (dash-dotted) 
is for the background-l(-2). Sep = +45° and the baseline is 
1,480 km. 



XII. DETECTOR SIZE AND GRANULARITY 

It is interesting to see what the effect of the detector 
size has on the performance of POL fit. Although the 
results reported so far are based on the analyses of the 
Monte Carlo events generated for the Super-Kamiokande 
detector with 40% PMT coverage, we can make some 
assessment of the effect of a larger, more granular de- 
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tector by imposing a cut on the distance to the PMT 
surface from the 7r° production point in the direction of 
the 7T° (D wa u). For this study we use single tt° events 
produced by NC interactions. When the ir° energy is 
about 1 GeV, the minimum opening angle of two pho- 
tons is about 20° and less at higher energies. As D wa ii 
gets larger, the number of PMTs that detect Chercnkov 
photons increases (improved granularity) which can help 
to resolve light patterns. Two negative impacts compete 
with this effect. The Cherenkov cone is more spread out 
so less light is seen by any one PMT. The light must 
travel further and is thus more subject to absorption 
and scattering (attenuation) in the water. These both 
lead to a decrease in the number of detected photons per 
PMT which may degrade the pattern of light that the 
Cherenkov cone produces. In the first case, since infor- 
mation is not lost and is just spread out to more PMTs 
the problem can likely be handled in improvements to re- 
construction codes. In the second case, some information 
about the original Cherenkov light emitting particles is 
unrecoverable. 

Figure [IT] shows the ir° detection efficiency as a func- 
tion of the opening angle for D wa u ranges: 5 m - 10 m, 
15 m - 20 m, and 25 m- 30 m and includes POLfit infor- 
mation. It is clearly seen that when the opening angle is 
smaller (less than 60°), the efficiency is improved as the 
distance to the PMT surface in the ir° direction increases. 
Note that above the opening angle of 60°, the ir° detec- 
tion efficiency seems more or less independent of D wa u. 
Because of the finite size of the detector and at such large 
opening angles selecting events based on D wa u becomes 
a less capable means of emulating a larger detector. 

Nevertheless, this indicates that the granularity of the 
detector in terms of the PMT density is an important 
factor to improve the 7r° detection efficiency. It also ap- 
pears that the effect of light attenuation is not a major 
issue at least for Cherenkov light traveling up to 30-40 m 
in SK-I water which has an attenuation length of 80-100 
m (depending on wavelength). 

Therefore for the same PMT coverage using the same 
PMTs, a larger detector than SK-I will perform better, 
on average, to reconstruct ir°. Limitations may be ex- 
pected once the typical path length for Cherenkov light 
approaches attenuation length. 

A similar study can be done by looking at how the 
S/B ratio varies as a function of D wa u for other event 
types. In this case we look at the S/B ratio for values 
of D wa u for the primary e-like ring. For D wa ii > 20 m 
the S/B ratio changes from the average of 1.4 to 3.8 for 
events with E rec < 1.2 GeV, and for events with 2 GeV 
< E rec < 4 GeV the S/B ratio essentially stays the same. 

The large improvement in the S/B ratio for events with 
Erec < 1-2 GeV results from an increase in the number of 
PMTs (pixels) in a Cherenkov ring. This improvement 
is significant as in the energy region E rec < 1.2 GeV the 
contribution from the NC events is reduced to a level 
as low as that from the irreducible background. This im- 
provement, however, is not realized for events with 2 GeV 
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FIG. 17: The n detection efficiency as a function of the two 
photon opening angle for three ranges of the distance from 
the 7T° production vertex to the closest PMT surface in the 
direction of n°. 

< E rec < 4 GeV presumably because in this energy re- 
gion multi-pion events are the major background. 

This observation is also true if the minimum distance 
D wa ii is set at 10 or 15 m, although the improvement in 
the S/B ratio is less than the case of the 20-m cut. For 
a SK-I sized detector, this 20-m cut reduces the number 
of the signal events by 41%. However, if the detector is 
larger, this loss of efficiency can be greatly reduced. In 
other words, for a given detector size, the finer granu- 
larity, not necessarily the number of Cherenkov photons 
collected by individual PMT, improves the S/B ratio. 

XIII. CONCLUSION 

For the baseline of 2,540 (1,480) km the detection ef- 
ficiency of the signal events using the SK-I cuts only is 
found to be 0.361±0.003 (0.373± 0.003) and the final ef- 
ficiency with the further 40% cut on the log-likelihood 
ratio is found to be 0.145±0.002 (0.149±0.002). The de- 
tection efficiency of the background- 1 events using the 
SK-I cuts only is found to be 0.054±0.001 (0.059±0.001) 
and the final efficiency with the further 40% cut on 
the log-likelihood ratio is found to be 0.0025±0.0002 
(0.0026±0.0003). The 40% cut on the log-likelihood 
should be considered as a guidance and the actual cut 
should be optimized after a detailed study of figure of 
merit that depends on the experimental design including 
the neutrino beam properties of a given experiment. 

We have demonstrated that a large water Cherenkov 
detector can be used to detect efficiently the signal events 
by v e s from the neutrino oscillation — > v e while keep- 
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ing the background at a reasonably low level in a VLBNO 
experiment for the baseline of over 1,400 km with a wide- 
band beam. 
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